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Neutrino oscillations in the early Universe prior to the epoch of primordial nucleosynthesis is 
studied by implementing real time non-equilibrium field theory methods. We focus on two flavors of 
Dirac neutrinos, however, the formulation is general. We obtain the equations of motion for neutrino 
wavepackets of either chirality and helicity in the plasma allowing for CP asymmetry. Contributions 
non-local in space-time to the self-energy dominate over the asymmetry for T > 3 — 5 MeV if the 
lepton and neutrino asymmetries are of the same order as the baryon asymmetry. We find a new 
contribution which cannot be interpreted as the usual effective potential. The mixing angles and 
dispersion relations in the medium depend on helicity. We find that resonant transitions are possible 
in the temperature range 10 < T <C 100 MeV. Near a resonance in the mixing angle, the oscillation 
time scale in the medium as compared to the vacuum is slowed-down substantially for small vacuum 
mixing angle. The time scale of oscillations speeds-up for off resonance high energy neutrinos for 
which the mixing angle becomes vanishingly small. The equations of motion reduce to the familiar 
oscillation formulae for negative helicity ultrarelativistic neutrinos, but include consistently both the 
mixing angle and the oscillation frequencies in the medium. These equations of motion also allow 
to study the dynamics of right handed as well as positive helicity neutrinos. 

PACS numbers: 13.15.+g,12.15.-y,11.10.Wx 



I. INTRODUCTION 



Neutrinos have taken center stage in particlephysics and have become a bridge between astrophysics, cosmology, 
particle physics and nuclear physics Lli Li LJ- LlJ- LJl • 

A wealth of experimental data have confirmed that neutrinos are massive and that different flavors of neutrinos can 
mix and oscillate 0, 0, S iiot - ITU Il2| thus providing indisputable evidence for new physics beyond the Standard 
Model. 

Neutrinos play a fundamental role in cosmology and astrophysics, and it is now widely accepted that resonant 
flavor oscillations due to the MSW effect in the sun provide a concrete explanation to the solar neutrino problem 

Neutrino oscillations in extreme conditions of temperature and density are an important aspect of Big Bang Nucle- 
osynthesis (BBN) and in the generation of the lepton asymmetry in the early Universe [3. IToL Il6l fill ITq . as well as in 
the physics of core collapse supernovae[rj ITH l2fj| . and the formation, evolution and cooling of neutron stars [2lll2^.l23| . 

While accelerator and reactor experiments measure the difference of the neutrino masses, high precision cosmological 
observations of the cosmic microwave background by WMAP combined with large scale structure suggest that the 
sum of the masses of all neutrino species is bound to be smaller than leV |24| . 

An important aspect of neutrino oscillations is lepton number violation, leading to the suggestion that the baryon 
asymmetry may actually originate in the lepton sector and the proposal that leptogenesis can be the main mechanism 
that explains the cosmological baryon asymmetry [2^, l2rl W5\ . 

Neutrino propagation in a cold medium has been first studied in ref. wherein the refractive index of electron 
neutrinos was computed. The early studies of neutrino propagation focused on the neutrino dispersion relations and 
damping rates in the temperature regime relevant for stellar evolution or big bang nucleosynthesis |l7l |28| . Since 
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then, the work has been extended to include leptons, neutrinos and nucleons in the medium [29|. The matter effects 
of neutrino oscillations in the early universe has been investigated in |30l l3l| . 

Since the original study of neutrino oscillations in the sunjlA IT^j. neutrino oscillations are typically studied within 
the single particle quantum mechanical formulation. For two flavors the evolution Hamiltonian is simply that for a 
two-state system, where the off-diagonal terms lead to the mixing and oscillation phenomena. The medium properties 
are input in this formulation after computing the contributions from charged and neutral currents in a medium 
with leptons, neutrinos and hadrons or quarks. A conceptually similar approach underlies the kinetic treatment of 
oscillations in the early Universe wherein the dynamics is studied from the time evolution of a density matrix that 
generalizes the single par ticle description but that does not generally account for the subtle aspects of flavor Fock 
states addressed in ref |37j] . which introduce a hierarchy of time scales. Recentlv|33| a set of generalized semiclassical 
Boltzmann equations for the single particle distribution functions, supposedly applicable to neutrino transport in core- 
collapse supernovae have been proposed, but where the mixing term must be obtained separately from the underlying 
field theory. For a thorough discussion of the kinetic approach to mixing and relaxation and the approximations 
involved, the reader is referred to [rH I32I l38l|. 

The dynamics of neutrino mixing in the presence of a background of neutrinos requires in gene ral a full non- 
linear treatment, which has so far been studied within a self-consistent single particle framework [34j in the form of 
approximate kinetic equations for the reduced density matrix [T^ l32| . Such study revealed a wealth of novel non- 
linear phenomena such as self-svnchronization|34j . An approximate treatment of background neutrinos within the 
framework of equilibrium finite temperature field theory has also been proposed in ref. [35]. However, a consistent 
treatment must necessarily rely on a non-linear kinetic description, which has not yet been developed in the full field 
theory framework. 

A calculation of the neutrino dispersion relations in a hot and dense medium implementing the techniques of 
quantum field theory at finite temperature was provided in ref. |28j| . This treatment was extended in ref. |35j to study 
the propagation of mixed neutrinos of negative helicity in a neutrino background, up to lowest order in g 2 /M^ 2 v . In 
ref.[2a| the quantum fields for Dirac neutrinos propagating in a cold but dense medium and the dispersion relation 
for both chirality components were obtained, again considering the medium corrections to the dispersion relations up 
to order g 2 /M^. 

The study of dynamics of neutrino mixing in the literature is mostly carried out within the framework of a single 
particle description, wherein the dynamical evolution is described in terms of an effective Hamiltonian for either a 
two or three level system (depending the number of flavors). However, a single particle formulation is inadequate in 
a hot and or dense medium where collective many body effects may be predominant. The main point in the above 
discussion is to highlight that there is a leap in the current approach to study neutrino oscillations in a medium: the 
result of a quantum field theory calculation of the index of refraction or effective potentials in the medium is input 
into a single particle quantum mechanical description of oscillations and mixing based on Bloch-type equations. 

We have previously reported on a study of oscillations in a hot and/or dense neutrino gas directly from the underlying 
quantum field theory[32j in free field theory. Such study, even at the free field level revealed quantum interference 
phenomena and subtle many body aspects responsible for a hierarchy of time scales that cannot be captured within 
the single particle description. In particular, the subtle aspects of the Fock representation of the distribution function 
as well as interference and coherence phenomena, leading to widely different time scales has not been fully included 
in the treatments in ref. 0, H2, H3 ■ 

More recently, novel collective neutrino excitations in the standard model (namely without masses and mixing) 
near the critical temperature were studied!^. Collective phenomena requires a systematic and consistent treatment 
implementing the methods of quantum field theory at finite temperature and density. Furthermore, the real time 
dynamics of mixing and oscillations requires a non-equilibrium formulation of quantum field theory specially suited 
to study the real time evolution as an initial value pro5/em|39(|-|43|. 

There are at least four fundamental reasons to study neutrinos in the early Universe at a deeper level: i) neutrino 
mixing may be the mechanism by which baryogenesis is a result of leptogenesi sl25l l26l |27| . ii) big bang nucleosynthesis 
is particularly sensitive to the spectrum and oscillations of neutrinosjlSt Il9l l30| . iii) just like the cosmic microwave 
background, there is a cosmic neutrino background left over from the big bang, and iv) neutrinos masses and mixing 
are the most clear experimental confirmation of physics beyond the standard model. All of these reasons warrant a 
complete quantum field theory study of neutrinos in the early Universe. 

The goals of this article: 

The full quantum field theory treatment of neutrino mixing in hot and or dense media has not yet received the 
same level of attention as the more familiar single particle treatment, which however, is not suited when collective 
many body phenomena become relevant as is t ypic ally the case in extreme environments. 

Previous quantum field theory studies|l7l l28l l35l l36| address either the dispersion relations or mixing phenomena 
under restrictive approximations to lowest order in g 2 /M^. 
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In this article we provide a systematic quantum field theory study of neutrino propagation and oscillations in the 
early Universe directly in real time. Because in the early Universe the lepton asymmetries are expected to be typically 
of the same order of the baryon asymmetry j?b/?7 7 ~ 10~ 9 a consistent description of neutrino propagation and 
oscillations requires to include corrections non-local in space-time of order g 2 /M w in the dispersion relations |28j and 
mixing angles. We focus our study on the case of two flavors of Dirac neutrinos, taken to be the electron and muon 
neutrinos, this study can be generalized to more flavors or to Majorana-Dirac mass matrices without any conceptual 
difficulty Our main goals are: i ) To provide a systematic and consistent study of the real time dynamics of neutrino 
oscillation and mixing directly in quantum field theory in conditions of temperature and lepton/neutrino asymmetries 
applicable to the early Universe prior to the nucleosynthesis era. This is achieved by formulating an initial value 
problem via linear response and implementing real time field theory methods at finite temperature and density, ii) To 
obtain the dispersion relations and in-medium mixing angles including the non local contributions from the neutrino 
self-energies up to order g 2 /Myy. The one-loop self-energy in expanded to lowest order in (oj, k)/Mw- We find a 
new contribution which cannot be interpreted as the usual effective potential. These contributions are necessary 
since the typical asymmetries in the early Universe are very small and these non-local (in space-time) contributions 
can be of the same order or larger than the local contributions. iii) To obtain the in-medium Dirac spinors for both 
helicities and study the evolution of oscillations and mixing for both helicity components directly in real time, iv) 
Two different temperature regimes are studied in detail: i) m e «T« m^, ii) m e ,m^ <C T <C Mw- The first regime 
is just prior to big bang nucleosynthesis. Lepton and hadron (proton and neutrons in nuclear statistical equilibrium) 
or quark asymmetries are included in the one-loop self-energy. We assess in detail the temperature and energy regime 
for which a resonance in the mixing angle is available in the medium. The second temperature regime is above the 
QCD phase transition and we include two flavor of (light quarks) with their respective asymmetries. In this regime 
the mixing angle becomes small. In both cases we also study the mixing and oscillations of positive helicity as well 
as right handed neutrinos, which are typically neglected in the literature. We also obtain the loop corrections to the 
oscillation frequencies thereby providing a complete description of oscillation and mixing that includes corrections 
to both the mixing angle and the oscillation frequencies, v) We obtain general oscillation formulae derived directly 
from the real time evolution in quantum field theory. These formulae reveal the limit in which the usual quantum 
mechanical single particle description is reliable as well as the corrections to them. 

Main approximations: Since our study relics on a one-loop self-energy computation including lcptons and neu- 
trinos, the inclusion of a neutrino background must necessarily imply some approximations for consistency. 

We do not yet consider absorptive contributions, which in the temperature regime studied here are of two- loop 
order, postponing the study of the interplay between oscillations and relaxation to a forthcoming article. 

Since we obtain the non-local (in space-time) contributions from the one-loop self energy we must address the issue 
of the neutrino propagators in the neutral current contributions. Because of mixing, the neutrino propagator in the 
flavor basis, in which the weak interactions are diagonal, does not correspond to the propagation of mass eigenstates 
and in principle the non-equilibrium propagators obtained in ref . |3?| must be used. The question of equilibration of a 
neutrino gas with mixing is one of time scales: the weak interactions are diagonal in the flavor basis, therefore weak 
processes tend to equilibrate flavor neutrinos with a typical weak interaction relaxation rate at high temperature 
[28j r ~ G 2 F T 5 . Oscillations, on the other hand mix flavors and tend to redistribute flavor neutrinos into mass 
eigenstates of energy E on a time scale r osc ~ E/5M 2 . Combined fitting of the solar and KamLAND data yield 
py| |<5M 2 | « 7.9 x 1CT 5 (eU) 2 , therefore considering E ~ T, we find T t osc ~ 10 (O.lT/MeV) 6 . This comparison of 
time scales suggests that for T > 10 MeV neutrinos are equilibrated as flavor eigenstates. Flavor eigenstates created 
at local weak interaction vertices will reach thermal equilibrium on time scales far shorter than those required for 
oscillations into mass eigenstates for temperatures larger than ~ 10 MeV. Since in a loop integral the typical momenta 
are of order T, and assuming the validity of this estimate, we consider the neutrino propagators in the neutral current 
self-energy loop to be diagonal in the flavor basis, massless and in thermal equilibrium. 

For temperatures T < 10 MeV and certainly below freeze out T < 1 MeV a full kinetic description that includes 
oscillations and expansion^, |3(| is required. The study of the kinetic equations will be the subject of forthcoming 
work. In this article we restrict our study to the temperature regime T > 10 MeV. 

We also assume that the lepton and neutrino asymmetries are of the same order as the baryon asymmetry, namely 
Li = (ni — n,i)/n 7 ~ 10~ 9 . For a relativistic species the asymmetry is proportional to & (1 + £j 2 /tt 2 ) with £j = fii/T, 
therefore under this assumption £j ~ 10~ 9 and we can safely neglect the contribution to the chemical potential in the 
non-local (in space-time) terms of order g 2 /M w . 

Brief summary of main results: The main results obtained in this article are: i) We obtain the equation 
of motion in real time for initially prepared wavepackets of neutrinos. Both chiralities and helicities are treated 
on equal footing. The self-energy is expanded up to lowest order in the frequency and momentum yielding non- 
local contributions to the equations of motion. These happen to be larger than those from the lepton and neutrino 
asymmetries for T > 5 — 10 MeV. ii) We studied two different temperature regimes: m e <C T <C m p within which 
we show that there is the possibility of resonant oscillations of test neutrinos, and m e , «T« Mw within which 
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the mixing angle for active neutrinos effectively vanishes. For T ~ fOMeV resonant flavor oscillations occur for 
neutrino energies in the few MeV range, iii) Mixing angles in the medium not only depend on energy but also on 
helicity. The dispersion relations of propagating neutrinos in the medium also depend on hclicity. Assuming that the 
lepton and quark asymmetries are of the same order as the baryon asymmetry in the early Universe, the non-local 
(in space-time) terms in the self-energies dominate over the asymmetry for typical energies of neutrinos in the plasma 
for T > 3 — 5 MeV. iv)The oscillation time scale in the medium is slowed- down near the resonance, becoming 
substantially longer than in the vacuum for small vacuum mixing angle. For high energy neutrinos off-resonance 
the mixing angle becomes vanishingly small and the oscillation time scale speeds-up as compared to the vacuum, v) 
The equations of motion reduce to the familiar oscillation formulae for ultrarelativistic negative helicity neutrinos, 
but they consistently include the mixing angles and the oscillation frequencies in the medium. We obtain general 
oscillation formula for cither chirality and helicity. These equations also describe the oscillation dynamics of right 
handed neutrinos, which, while suppressed as consequence of the small masses, are not sterile. 

The article is organized as follows: in section [n] we obtain the equations of motion for initially prepared neutrino 
wavepackets by implementing the methods of non-equilibrium field theory and linear response. In section IIIII we 
obtain the one-loop self-energy contributions from charged and neutral currents. Section llVI is devoted to obtaining 
the dispersion relations, mixing angles and oscillation time scales in the medium and a study of the possibility of 
resonances. In section^we study the real time evolution of neutrino wavepackets as an initial value problem. Section 
I VII presents our conclusions, summarizes our results and presents some conjectures and further questions. The detailed 
calculation of the self-energy is presented in an appendix. 

II. EFFECTIVE DIRAC EQUATION FOR NEUTRINO PROPAGATION IN A MEDIUM 

The propagation of a neutrino in a medium is determined by the effective Dirac equation which includes the self- 
energy corrections. Its solution yields the real time evolution as an initial value problem. The correct framework to 
study the dynamics is the real time formulation of field theory in terms of the closed-time-path integral [3^ -^3- ^ n 
this section we implement this method combined with linear response to obtain the effective equation of motion for 
an expectation value of the neutrino field. The main concept in this approach is the following, consider coupling an 
external c- number Grassman source to the neutrino field and switching this source adiabatically up to time t = 0. 
This source induces an expectation value of the neutrino field, after switching-off the external source at t = 0, the 
expectation value evolves in time as a solution of the effective Dirac equation in the medium with the initial condition 
determined by the source term. 

The main ingredient in this program is the retarded self-energy which enters in the effective Dirac equation. The 
real-time formulation of field theory directly leads to causal and retarded equations of motion. It is important to 
highlight the difference with the S-matrix approach which describes transition amplitudes from in to out states, the 
real time formulation yields the equations of motion for an expectation value and these are fully causal 39]- 43]. 

The self-energy is obtained in the unitary gauge in which only the correct physical degrees of freedom contribute 
and is manifestly unitarvflil]. Previous calculations of the neutrino self-energy in covariant gauges (one of which is 
the unitary gauge) have proven that although the self-energy does depend on the gauge parameter, the dispersion 



As mentioned above we restrict our discussion to the case of two flavors of Dirac neutrinos, namely the electron 
and muon neutrinos. The subtle CP violating phases associated with the case of three active neutrinos will not be 
considered here. However, the method can be generalized to three active neutrinos, sterile neutrinos or even Majorana 
neutrinos without any conceptual difficulty and will be postponed for further discussion elsewhere. 

For Dirac neutrinos, mixing and oscillations can be implemented by a minimal extension of the standard model 
adding a Dirac mass matrix to the standard model Lagrangian which is off-diagonal in the flavor basis. The relevant 
part of the Lagrangian density is given by 




C — C v + C w + C z 



+ C-CC + £-NCi 



(III) 



where £ 



is the free field neutrino Lagrangian minimally modified to include a Dirac mass matrix 



C v = v a (i$5 ab 



M ab ) v b 



{11.2) 



with a, b being the flavor indexes. For two flavors of Dirac neutrinos the mass matrix M ab is given by 




(II.3) 
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C^z are the ^ Tee fi e ld lagrangian densities for the vector bosons in the unitary gauge, namely 

c° w = -~ (^w+ - 9,w+) " - d v w- ") + w+ w~ " , (11.4) 

4 = -J (d^ - c^) - d v Z») + ijlff Z p Z» , (II.5) 

and the charged and neutral current interaction lagrangian densities are given by 

CcC = ^= [Va r L l a W+ + Z o7 " L V a W~] , (II.6) 

Cnc = TT^fT [ Fq ^ L Va Z » + 7* ^ (Sa ~ 9t 7 5 ) fa Z„] ■ (H.7) 

Zi COS U yj 

where L = (1 — 7 5 )/2 is the left-handed chiral projection operator, g V ' A are the vector and axial vector couplings for 
quarks and leptons, I stands for leptons and / generically for the fermion species with neutral current interactions. 

For two flavors, the diagonalization of the free field Dirac Lagrangian for neutrinos, 1II.2|) is achieved by a unitary 
transformation to mass eigenstates. Considering, flavor and mass doublets respectively 



related by unitary transformation 



with the unitary transformation given by the 2x2 matrix 



= tff !? ) . ax 



sm t/ cos f 

where is the vacuum mixing angle. 

In the basis of mass eigenstates (y\, 1/2) the mass matrix M ao becomes diagonal 

Mi 
M 2 



The elements m^m^ and m eAt in the mass matrix (|II.3|I are related to the vacuum mixing angle 9 and masses of 
the propagating mass eigenstates Mi and M2 as follows 

m ee = C 2 Mi + S 2 M 2 ■ m w = S 2 M x + C 2 M 2 ; m ejU = -(Mi - M 2 ) C S , (11.10) 

where C = cos 8 and 5 = sin 0. 

For later convenience, we introduce 

M=Ml±Ml . 5M 2 =M?-Mt (11.11) 

The current value for the average of the vacuum masses obtained by WMAP|24| and oscillation parameters from the 
combined fitting of the solar and KamLAND data are |45j : 

M«0.25(eF) ; |<5M 2 | w 7.9 x 10 -5 (eV) 2 ; tan 2 9 « 0.40 . (11.12) 

For these values of the masses and more generally if there is an almost degeneracy in the hierarchy of neutrino masses 
the ratio 

J^l « 1 . (...13) 

M 

The smallness of this ratio in the nearly degenerate case will lead to important simplifications. 






FIG. 1: One-loop diagrams contributing to the neutrino self-energy. 



Our goal is to obtain the effective Dirac equation for neutrinos propagating in the medium and extract the in- 
medium mixing angles, propagation frequencies and the wave functions of the propagating modes in the medium. 
The real-time effective Dirac equation in the medium is derived from linear response by implementing the methods 
of non-equilibrium quantum field theory described in |43| . 

Following this approach, we introduce an external Grassmann-valued source that couples linearly to the neutrino 
field via the lagrangian density 

Cs = v a r) a + r] a v a . (H.14) 

whence the total lagrangian density is given by C + Cs ■ The external source induces an expectation value for the 
neutrino field which, in turn, obeys the effective equation of motion with self-energy modifications from the thermal 
medium |43j |. 

To study the dynamics of the system, it is the expectation values rather than the in-out S-matrix elements that are 
necessary 40J. This requires a generating function for the real-time correlation functions. Denoting generically by $ 
the fields (fermions, or gauge bosons), a path integral representation of this generating functional is given by 

Z\j+,r] = J DQ+DQ-eWP+^-W*^ , (11.15) 

and the path integrations over the fields < I> ± will be taken along the forward (+) and backward (— ) time branches, 
in the presence of the sources [3|j-^i|. Here, the sources are coupled linearly to the fields $ ± and thus 
the real-time correlation functions can be obtained from the functional derivatives of this generating functional with 
respect to these sources. Functional derivatives with respect to j + and j~ give the time-ordered and anti-time-ordered 
correlation functions respectively. 

The sources are introduced to compute the real-time correlation functions and will be set to zero after the 
calculations. However, the external Grassman source r\ and hence the expectation value induced to the neutrino field 
will remain the same along both time branches (±). For further discussions on the general method, see references 

The equation of motion for the expectation value of the neutrino field induced by the external Grassman source is 
derived by shifting the field 

vt=i>a + ^t , 1>a = (v£) , <^)=0, (11.16) 



and imposing (Vf^) = order by order in the perturbation theory |4ll l42t l43|. Carrying out this implementation up 
to one-loop order, we find the following equation of motion 



(ipSab-Mab + T,^ L) Mx,t) + J d 3 x' J dt'Z r a ?(x-x',t-t')MZ',t') = -ria(Z,t), (11.17) 

where E£|*(af — x', t — t') is the real-time retarded self-energy given by the exchange one-loop contributions, the first 
two diagrams displayed in Fig. H an d L is the tadpole contribution from the last two diagrams in Fig. ^ The 
expectation value of the neutrino field in the medium describes a beam or wave-packet of test neutrinos, namely these 
are neutrinos that are injected in the medium, for example by the decay of a neutron or any other heavy particle, but 
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have not (yet) thermalized with the plasma. This is precisely the manner in which linear response leads to a study of 
real-time phenomena. 

Due to the translational invariance of the medium in thermal equilibrium, the retarded self-energy is simply a 
function of x — x' and t — t'. Hence, the equation can be written in frequency and momentum space by introducing 
the space-time Fourier transform of the expectation value and the source 

Mx,t) = ^=J2j duj Mk,u) e iS - s e"^ 4 , Va (x,t) = J= ^ J doj Va (k,uj) (11.18) 

k k 

Furthermore, due to the rotational invariance of the thermal medium, implies that all physical quantities depend on 
\k\ = k. 

We have argued in the introduction that for temperatures T > 10 MeV the relaxation via the weak interaction is 
faster than the time scales of (vacuum) oscillations. The validation of this assumption requires a deeper study of 
the interplay between neutrino oscillations and damping processes in the medium including the two-loop self-energy, 
which will be the subject of a forthcoming article. Assuming the validity of this estimate, the neutrinos in the loop 
are in thermal equilibrium as flavor eigenstates and at the temperatures of interest we consider them massless. Under 
this assumption, the self-energy is diagonal in the flavor basis 

As a result, the effective Dirac equation for neutrino oscillations in the medium is obtained as 

( 7 ° u - 7 • k) S ab - M ab + Tft* L + E Qb (w, k) L] Vb(w, *) = k) , (11.19) 

with 

E(w,fc) = + Y, z {to,k) ; E tad = T, Ve tad + E w " tad + E-^ tad , (11.20) 

where Y,f tad is the total tadpole contribution from all fermions other than neutrinos in the loop. 
The details of the calculation of the self-energy are presented in the Appendix. 

III. ONE-LOOP SELF-ENERGIES 

A. Charged and Neutral Currents 

From Eas. (|A.28(l - (|A.32() in the Appendix, we found that the self-energies k) can be written in the dispersive 

form 

f dk ImE w (k ,k) 
Y lW ,z{u,k)= I — . (HIT) 

J 7r ko — u> — le 

with 

/mEuKfco, k) = J | [1 - Nf(u p ) + N b (W q )} fQ(p, q) S(k -lo p - W q ) + 

+ [1 - N f (u p ) + N b (W q )] @(-p, -q) 5(k Q +uj p + W q ) + [N f (u p ) + N b (W q )} Q(p, -q) 5(k - u p + W q ) + 

+ [Nf(u p ) + N b (W q )] @{-p, q) 5(k + uj p - W q ) 1 , (III.2) 



M w 



where we have defined 



qK = (W qi k-p) , to p = ^\p\i+m} , W q = ^\k-f\^+M^. 
The corresponding contribution from neutral currents can be obtained from the above expression by setting 

• M W ^M Z = ^- , (IH.4) 
sjl 2 cos W cos W 
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with 9 W being the Weinberg angle. 

In the limit T -C Myy z-i the abundance of vector bosons is exponentially suppressed, hence we neglect the terms 
that feature N (W q ). The imaginary part of the one- loop self-energy vanishes on the neutrino mass shell at one-loop 
level. A non-vanishing damping rate (non-vanishing imaginary part of the self-energy at the neutrino mass shell) at 
temperatures T -C Mw arises at two-loop level. Thus we focus solely on ReY>w(ko, k) when studying the dispersion 
relation and propagation of neutrinos in the medium. 

The form of ImT,w(ko, k) suggests that ReY>w{^, k) can be written as 



Re^ w (io, k) = 7° OyyiU) k) - 7 ' k ^wi^, k), 



(III.5) 



where ct^(cj, k) and a^ v (w, k) can be obtained by taking traces on both sides. 

Dropping the T = part and using the dispersive representation wc find that for any fermion / in the loop, 



d 3 p 



2 J (2tt) 3 4 W q io p 




Qo(p,q) Qo(p,-q) 



W q + UJp — U) W q — UJ p + Uj) 

Qo(p,q) Qo(p,-q) 



W q + Ldp + UJ W q — UJp — UJ 



(III.6) 



d 3 p 



v w ^k) = --J — 



) 3 4 Wq Up 



N f (u p ) 



k • Q(p, q) k Q(p, -q) 



W„ + Ld v - UJ W 



+Nf(u P ) 



k-Q(p,g) k-g(p, -q) 



W n 



w n 



(III.7) 



In the thermalized medium with temperature T, the dominant loop momenta is of order p ~ T, therefore we neglect 
the neutrino masses since T > 3 — 5 MeV. The self-energy is expanded in a power series in the external frequency and 
momentum ui, k, we refer to terms of the form uj/M-w ; k/Myy as non-local terms (in space-time) since they represent 
gradient expansions in configuration space. Furthermore, we neglect the contribution from leptons with masses rrif 3> 
T since these will be exponentially suppressed, but we calculate the self-energies up the order (g 2 jM^ z )(rrif /T) 2 for 
leptons with masses m/ -C T and all higher order terms will be dropped. A straightforward but lengthy calculation 
gives 



r (u),k) = - 



>W 



3G F 
7tt 2 



G F 



n 



f 

k T 4 



) + 



7tt 2 G f uT 4 



45 y/2 



Ml 



15 V2 
30 



1 - 



/m/V 
7 7T 2 V T J 



(LIU 



where Gf = V% 
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any fermion / defined as 



is the Fermi constant, and nf — tit is the particle-antiparticle number density difference for 



= 2 



d 3 p 
(2tt) 3 



(III.9) 



The contribution to crh^ioj, k) of order g 2 /M^, vanishes as a consequence of the isotropy of the equilibrium distribution 
functions. In calculating the non-local (in space-time) terms proportional to lu/Mw ; kjMyj we have neglected the 
chemical potentials under the assumption that all asymmetries (leptons and neutrinos) are of the same order as the 
baryon asymmetry, in which case fi/T ~ 10~ 9 where [i is the chemical potential for the corresponding species. 

In ref.|28j the equivalent to er° is quoted as the coefficient 6l, but no equivalent to a 1 was provided there, this is a 
difference between our results and those of ref.|28l|. The contribution u 1 ^, k) is new and it cannot be identified with 
an "effective potential" (which is proportional to 7 ) and it will be a source of helicity dependence on frequencies and 
mixing angles which has not been appreciated in the literature. 

To obtain the corresponding expressions for the neutral current interactions, we can simply apply the replacement 
The tadpole contributions have been previously obtained in ref. 28] where we refer the reader for more 



ca. (|lTL4| 

details. The result of the tadpole diagrams is summarized by the expression 



tad 



2 Gf g 
V2 



v 



(nf - nj) 7 , 



(111.10) 



where the coefficient g v for the different species of leptons, hadrons and quarks are given in the table below. 
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Particles 


9 V 


Ve, U/i, V T 


1 

2 


e,n,T 


1 
2 


h 2 sin 2 0„, 


P 


1 
2 


2 sin 2 0„, 


n 


1 
2 


u, c, t 


1 

2 


| sin 2 e w 


d, s, b 




- | sin 2 6»„ 



Writing Z ftad (uj, k) in the same form as (|III.5|) . we obtain 

o 2G F g v 



'tad 



V2 



(n/ - n/) 



-'tad 



(III.ll) 



and tadpoles with quark loops acquire an extra factor 3 from color. The tadpole contribution is proportional to 
rif — rit, which is the signature of a CP asymmetric medium. It is customary and convenient to express rif — nj in 
terms of its relative abundance to the photons in the universe. At any temperature T, the photon number density in 
the universe is given by 



n, = - C (3) T 3 . 

and the particle-antiparticle asymmetry for any fermion species / is defined as 



(III. 12) 



(111.13) 



The magnitude of the observed baryon asymmetry is Lb — 10 -9 . Since B — L is conserved in the standard model, 
it is natural to expect that the lepton asymmetries should be of the same order as Lb- Although there is no a priori 
reason to expect the neutrino asymmetries to be of the same order, we will henceforth assume that all the lepton and 
neutrino asymmetries L ei L^, L Vp and L v are of the order of the baryon asymmetry ~ 10 -9 . 

The next step is to compute the contributions from the one- loop exchange diagrams. We focus on two different 



temperature regimes, m e «T« uij, in which only the electron is ultrarelativistic and m e 
both leptons are ultrarelativistic. 



«T« Mw in which 



1, m e < T < m M 

This temperature limit is interesting because it is the energy scale right above BBN. The efficiency of BBN is 
sensitive to the amount of electron neutrinos which, in turn, depends on the detailed dynamics of neutrino oscillations. 
As it will be discussed below in detail, in the temperature regime T > 5 McV the non-local (in space-time) terms 
proportional to uj, k from the exchange diagrams (both charged and neutral currents) dominate the contributions of 
the lepton and neutrino asymmetries assuming all of them to be of order 10 -9 . 

In the temperature limit with m e <T-C m^, the contribution from fi leptons to the exchange one loop diagram is 
exponentially suppressed and we neglect it. Apart from electrons and neutrinos, the thermal background does contain 
protons and neutrons in nuclear statistical equilibrium since for T > 1 MeV the weak interactions lead to equilibration 
on time scales shorter than the Hubble time scale via the reactions 



n <-» p + e 



p + e <-> n + v e 



; e , p + v e «-> n + e T 
In the basis of flavor eigenstates, the total one-loop self-energy contribution is of the form 

Re k) = [7°A(w) - 7 • kB(fc)] L 
where A(w, k) and B(w, k) are 2x2 diagonal matrices in the neutrino flavor basis given by 



A(u) 




»(*) = 




(III. 14) 



(111.15) 



(111.16) 



10 



Extracting non-local terms (in space-time) up to G(ui/T) , 0(k/T) we find the following matrix elements, 

.2 n-,,rr4 qn„ 7 ^2 g f ojT 4 



A e (ui) 



3 G F 7 7r Gf w T 

G F 



15 M? 



w 



4 3G F 

2\/2 30 V2 



V2 



—j= L Vfl n 7 + ——(-- + 2 sin 6 W ) L e n 7 



2 G F / 1 • 2 A \ t F 



(III. 17) 



3 Gf , 

+ 7f * 711 + 

2 G P ,1 



7 7T 2 G F w T 4 



30 V2 



Ml 



2 C 1 

^_ n 7 + — + 2 sin 2 W ) L e n 7 



C, 



H 7S~(o - 2 sin L p n 7 -= L„ 



B e {k) 



V2 V 2 
7 7T 2 G F k T 4 



45 V2 M 2 , 
7tt 2 G F kT 4 
' 90 V2 M 2 



1 - 



V2 

30 /m e ^ 2 ' 
T 



'7 ' 



7tt 



7 tt 2 G F fc T 4 



90 V2 



M 2 



(111.18) 
(111.19) 
(111.20) 



We purposely displayed the individual terms in the above expressions to highlight that the first line in A e ^{uj) as 
well as the expressions for B e ^(k) arise from the exchange diagrams (the two top diagrams in figQ, while the second 
and third lines in A e ^{uj) arise from the tadpole diagrams (bottom two diagrams in fig^). We have assumed that 
the flavor neutrinos are in thermal equilibrium and have consistently neglected neutrino masses. The correction term 
(m e /T) 2 is displayed so that one can estimate the error incurred when this term is dropped for T » m e , the error 
is less than 1% for T > 5 MeV. In what follows we will neglect this contribution in the temperature range of interest 
for this section. 

Charge neutrality requires that L e = L p , hence 



tad 



Y^p tad q 



Therefore, the expressions for A e and A^ simplify to the following: 



A e (u) 
B e (k) 



Gf n 7 
Gf n 7 



Ce 



7tt 4 uoT 



60C(3) Af 2 , 
7tt 4 ujT 



'2 + cos 2 0„ 



60C(3) M 2 , 



cos 



Gf fin, 7 7r 4 k T , n 

— —f^ 1 r~r — o- (2 + cos 2 1 

V2 180 C(3) M 2 , v 



Gf 7 7r 4 k T 
V2 180 C(3) 



cos 2 9 W . 



(111.21) 

(111.22) 

(111.23) 

(111.24) 
(111.25) 



w 



where 



- L Vtl — 3 L e — L n 
L v — L n 



(111.26) 



Where in the temperature regime m e <C T <C we consistently neglected the muon contribution to the non-local 
terms proportional to uj, k in the exchange diagrams. 

The above expressions also reveal the importance of the temperature region T > 5 MeV. Assuming that all 
asymmetries are of the same order as the baryon asymmetry, namely Li ~ 10 -9 , we see that for u> ~ k ~ T the factor 
T 2 /M^ > 10~ 9 for T > 5 MeV. 

We will discuss below that in this region there is also a resonance in the mixing angle in agreement with the results 
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2. m e ,m M <T<M w 

This temperature region is important because the non-local contributions arc much larger than that of the lcpton 
and neutrino asymmetries, assuming both to be of the same order ~ 10~ 9 and are the same for both leptons if their 
masses are neglected. Therefore, if the contribution from the lepton and neutrino asymmetries is neglected, and 
terms of 0(m 2 /T 2 ) ; 0(m^/T 2 ) are neglected, the matrices A, B become proportional to the identity. In this case the 
mixing angle would be the same as in the vacuum. We will see however, that keeping terms of 0{m 2 e /T 2 ) ; C(m^/T 2 ) 
leads to a very different result, namely the vanishing of the mixing angle for negative helicity neutrinos or positive 
hclicity antineutrinos in this temperature range. 



For T 3> m b 



100 MeV the temperature is larger than the critical temperature for deconfinement in QCD 



T c ~ 160MeV. Therefore, the medium contains free quarks but no nucleons. Since both u and d quarks have masses 
smaller than 10 MeV, their masses can be neglected. We only include in our description the two lightest quark degrees 
of freedom consistently with keeping only a weak doublet. We also assume that there is vanishing strangeness in the 
medium. 

As a result, the corresponding A e:tl (u), B e ^{k) are now given by 



Ae(u) 



3 Gf 7tt 2 G f luT 4 

G F 



15 y/2 M 
G F t 

+ "7f Uj + ni 



3 G F 
2y/2 



7 7T 2 G F LO T 4 



30 y/2 



M 2 



+ 



+ 



2 Gf , 1 „ . 9 „ % T 2 Gf , 1 ■> 

—^-(-- + 2 sin 2 8 w )L e n 1 + ^' 



(-| + 2 sin 2 9 W ) n 7 



6G f .l 4 . - . T 

(o ~ o Sm °™> L u n J + 



V2 v 2 
3 Gf t 



7 ir 2 G F uT A 



15 V2 M 2 



6 Gf 1 1 2 . 2/1 . _ 

~^y~2 + 3 sm ' 711 ' 

3G F 7 7r 2 Gf w T 4 



(111.27) 



w 



2 V2 



30 



+ GF L 



I ^ F T 

n 1 + — U 



2 1 o 2 Gp . 1 

+ ^7T ( ~2 + 2 w) Le " 7 + ~V2~ ( ~2 + 2 sin 

6G F/ 1 4 . 2n . T 6G F/ 1 2 

+ ^7T ( 2 - 3 Sm M iu % + "7T ( 2 + 3 Sn: 



Be(*0 



7 7T 2 Gf k T 4 
45 v/2 M 2 , 

7 7T 2 G F fc T 4 
45 V2 MA, 



1 



1 



30 /m e y 

7 7r 2 v r / 



) « 7 

; ' 3 ) id % , 
7 7T 2 Gf fc T 4 



30 
7^ 



T 



90 >/2 M 2 
7 7T 2 Gf k T 4 
90 V2 Aff 



(111.28) 
(111.29) 
(111.30) 



Charge neutrality of the medium leads to the constraint 4 L u — Ld — 3 L e = 0, which leads to the following simplified 
expressions: 



Ae(w) 
Be (A) 



Gf n 7 
Gf Tij 
Gf n 7 



"£ e + 7^7X7 TJT I 2 + COS 2 0„ 

M w 



C 

7 ir 4 



60 C(3) 
7tt 4 ujT 



" 60 C(3) M 2 , 



cos 



V2 180 C(3) M w L 
Gf n 7 7 7r 4 



fcT 



V2 180C(3)Af, 2 



2 + cos 2 e v 
2 + cos 2 6 V 



60 
7^ 
60 



i^ 2 \~) 



(111.31) 
(111.32) 

(111.33) 
(111.34) 



where 



-Ce = --L, e +^-3I e + (l-4 sin 2 9 W ){2 L e - L M ) - (1 - 8 sin 2 6> W ) J L„ - 2 L d 



L Ve - 3 L„ + (1 - 4 sin 2 fl w )(2 L e - L„) - (1 - 8 sin 2 6 W )L U -2L d 



(111.35) 
(111.36) 
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In the limit when T 3> m e ^ both leptons become ultrarelativistic and a CP-symmetric medium becomes flavor blind 
to the weak interactions. In this case we must keep terms of 0(m etfM /T) to understand the nature of oscillations and 
mixing. 



IV. DISPERSION RELATIONS, MIXING ANGLES AND RESONANCES IN THE MEDIUM 

The neutrino dispersion relations and mixing angles in the medium are obtained by diagonalizing the homogeneous 
effective Dirac equation in the medium, namely by setting r)(u>, k) = in ca. (|II.19(l . Using the results obtained above, 
the homogeneous Dirac equation in frequency and momentum becomes 



7° tut - 7 • k kt - M + (7 A(w) - 7 • kB(k)) L ip(u, k) = , 
where 1 is the 2x2 identity matrix in the flavor basis in which the field ip(u>,k) is given by 

v e (uj,k) \ 



tp(w,k) = 



t {u),k) 



(IV.l) 



(IV.2) 



with v e (u>, k) and ^(w, k) each being a 4-component Dirac spinor. 

If we multiply the effective Dirac equation l|IV.l|) by the chiral projectors R and L respectively from the left, we 
obtain 



(7 W-7- kl) ip L -Mip R = , 
(7 u> 1 - 7 • k k tj ip R - M ip L = 



where we have defined the flavor matrices 



K = kt 



(IV.3) 
(IV.4) 

(IV.5) 



The set of equations (IIV.3|) and i|iV4|) couple ipL and ipR together. To solve the equations, we first multiply l|IV.3|) 
by ^7° uj — 7 • \tkj 1 from the left and use ea. (jIV.4|l to obtain an equation for ipL which can be written in terms of 
the helicity operator /i(k) = 7°7 • k7 5 , as follows 

wW-H + ii(kp-Wi)-M 2 V>l = 0, 

and the right handed component is given by 



(IV.6) 



^,fc)-M 7 ° [ " + 2 M1 : ) 2 fc] M",k) 

lo z — k z 



(IV.7) 



It is convenient to separate the Dirac and flavor structure to simplify the study. This is achieved most economically 
in the chiral representation of the Dirac matrices, in which 



-1 

-1 



7 5 = 



1 

-1 



(IV. 



7 • k = 



akt 



h(k) = a- k 



1 



v -<r-kl / ' v ' V 01 / 
and by introducing the two component Weyl spinors w^^(k) eigenstates of helicity, 

a-kv ih] (k) = hv ( - h) (k) ■ h = ±l. 

These spinors are normalized so that 

(«W(£)) t vM(k) = 6 h , h , . 



(IV.9) 



(IV.10) 



(IV.ll) 
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In terms of these helicity eigenstates, a general flavor doublet of left (L) and right (R) handed Dirac spinors can be 
written 



/l — ±1 \ 



and 



where ip^ ; £W are flavor doublets. We have purposely left the arguments unspecified because this expansion will be 
used in real time as well as for the Fourier and Laplace transforms respectively. We need both positive and negative 
helicity eigenstates because the four independent degrees of freedom for each flavor are positive and negative energy 
and positive and negative helicity. 

Projecting eq. ljIV.61) onto the helicity eigenstates ^^(k) we obtain an equation for the flavor doublet tp^(u>, k), 
namely 

[{uj 2 - k 2 )t + (w - hk)(A + MB) - M 2 ] if {h \uo, k) = . (IV.14) 
Projecting eq. (|IV.7() onto helicity eigenstates yields the relation 

= — M (tJ 2 +/t | ? ¥> (& W) . (IV.15) 
uj z — k l 



Writing the doublet fc) in the flavor basis as 



^k)=[^l 



a™ b \ ( v&> 
b C W ) \ v f> 



(IV.16) 



leads to the following matrix form for Ea. (|IV.14|) 

= , (IV. 17) 



where the matrix elements in the flavor basis are given by 

a (h) = uj 2 - k 2 + (ui - hk)(A e + hB e ) - ~ {M\ + Aff ) - ~ (M\ - M%) cos(26>) , (IV.18) 

b = - (Aff - M 2 2 )sin(2(9) , (IV.19) 

c {h) = uj 2 -k 2 + (w - hk){A„ + hB^) - i (A/ 2 + Af|) + ~ (A// 2 - Aff) cos(26>) . (IV.20) 
Let us introduce a doublet of mass eigenstates in the medium 

x^,k)=( v f^\ f (rv . 21) 

related to the flavor doublet ip^ {uj, k) by a unitary transformation uj$ with 

CTW = f C0S ^ } S[ne Tl) , (IV.22) 
\ - sin 9 m cos 6m J 

<pl h \u,k) = U<Px (h Hu,k) ; (W(^fc) = f/W(W( W) i), (IV.23) 



14 



The mixing angle in the medium for states with helicity h, 9$ is obtained by requiring that the unitary transformation 
cq. (|IV.22|) diagonalizes the matrix equation ea. (|IV.17|) . The eigenvalue equation in diagonal form is given by 



J 2 - k 2 + \ S h (uj, k) - I {Ml + M 2 ) -±6M 2 



where Sh(u>, k), SM 2 and are respectively given by 



\26 



A h (u>,k) 
SM 2 



sin 2 29 



1 

-1 



S h (uj,k) = (oj -hk)[A e (uj) + Af,{uj) + h B e (k) + h B^k)} 

SM 2 = M'l - Ml , 
A h (w,k) = {uj-hk)[A e {uj)- A^{uj) + hB e {k)-hB^k)] 

The mixing angle in the medium is determined by the relation 

2 b SM 2 sin(26») 



tan[2flW] 



z (h) _ a (h) SM 2 cos (26») - A h (u, k) ' 



or alternatively by the more familiar relation 

sin 26^ = 



sin 29 



COS/V jjp- 



sin 2 29 



(IV.24) 

(IV.25) 
(IV.26) 
(IV.27) 



(IV.28) 



(IV.29) 



We note that the neutrino mass eigenvalues as well as the mixing angle depends on k as well as on the helicity 
eigenvalue h. This is one of the novel results which has not been obtained before simply because only left handed 
negative helicity neutrinos were considered in the literature P, ITTl [T3. fl3L fill Il5| . 

The right handed components are obtained from the left handed ones by performing the unitary transformation 
eq. on eq. (|IV.15|) . The relation HIV. 15)1 leads to the following expressions 



k 2 



SM 2 C S 



AM 



2 1 s -c 



X {h) (oj,k) 



where M = h(Mi + M 2 ) and 



(IV.30) 



C = cos 



2 0#) - 2 



S = sin 2 - 2 



(IV.31) 



A. Eigenvectors and dispersion relations 
Ea. l)IV.24p has the following eigenvectors in the basis of mass eigenstates: 



X[ h \u,k) = V [ n \u,k) 



and 



Cf>>*0 



.;)■ 

lu + hk — 
M 



j 2 -k 2 



-v\ (Ld,k) — — M 



J 2 -k 2 



SM 2 I C 

aW [ s 



SM 2 
AM 2 



S 

-C 



(IV.32) 
(IV.33) 

(IV.34) 
(IV.35) 



The corresponding doublets in the flavor basis can be obtained by the unitary transformation ea. l)IV.22)) . 
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The eigenvalues are found in perturbation theory consistently up to 0(Gf) by writing 



, a=l,2 



with 
We find, 



E U2 (k) = Jk 2 + Ml 



m(k) 



S h (±E 1 {k),k)~5M 2 



--Mil 



- 1 



(IV.36) 



(IV.37) 



(IV.38) 



6J 2 h \k,±) = - 



4£ 2 (fc) 



S h (±E 2 (k),k) + SM 2 



cos 29 



A h (±E 2 (k),k) 
SAP 



sin 2 26> 



- 1 



(IV.39) 



It is important to highlight that whereas the mixing angle only depends on Ah, the medium corrections to the 
frequencies also depend on Sh- This is important because even when in the case when the matrices A, B become 
proportional to the identity, in which case = and the mixing angle in the medium coincides with that of the 
vacuum, the frequencies and in particular the oscillation frequency still receives medium corrections. 



B. Resonances 

The condition for resonant oscillations is that the mixing tan[2^m' ) ] reaches a maximum (infinity) as a function of 
a parameter, temperature, density or energy. From eq. i|lV.28[) a resonance takes place when 

0V.40) 

where to — ui h \k,±) correspond to the dispersion relations for the propagating modes in the medium, given by 
cq. (|IV.36|) . To leading order in Gf the in-medium dispersion relation can be approximated by the free field dispersion 

relation lJ£^ (k, ±) « ±^Jk 2 + M 2 2 . The relativistic limit is warranted because the neutrino momenta in the plasma 

is k ^S> M\ t 2 ~ eV. Furthermore, under the assumption that the hierarchy of vacuum mass eigenstates is nearly 

degenerate, namely \SM 2 /M | <§; 1, as seems to be supported by the experimental data, the dispersion relations can 
be further approximated as follows 

W W(fc,±)«A*^l + ^ ; A = ±l (IV.41) 

It is convenient to introduce the following notation 

C 9 - lO 9 ^-^) (IV.42) 



eV" 



S* = 105 ( ) (IV-43) 



// the lepton and neutrino asymmetries are of the same order of the baryon asymmetry, then 0.2 < |£g| < 0.7 and 
the fitting from solar and KamLAND data suggests \S$\ w 8. Using the approximations leading to eq. (|IV.41fl the 
ratio A h /5M 2 can be written compactly from ea. (|lV.27|) . 

We study separately the cases m e <C T and m e , m M <C T <C Mw . 
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1. m e «T« m,, 



case I: ui = k + 



M 
2k 



-1, positive energy, negative helicity neutrinos: 



A h _ £ /0.1TV 
<5M 2 ~ <5 5 VMeV / T 



-Co 



( 2T 



k 



VMeV J T 



(IV.44) 



Where we have neglected =p-. For fixed temperature, the resonance condition ea. (|IV.40(l is fulfilled for the value 
of neutrino momentum given by 



(MeV) 2 J 



^ + 16. 5 cos2,(^X) ; 



(IV.45) 



Hence, for 5$ cos 20 > 0, there is always a resonance. If \Cg\ < 1, then for neutrino momenta such 
that VTk > 1 MeV the non-local term dominates over the asymmetry and the resonance occurs for k ~ 
25 v 7 ^ cos 29 (MeV) 3 /T 2 . For example, if T ~ 10 MeV, the resonance occurs for k ~ 1 MeV. If 5 5 cos26> < 
there can also be a resonance provided 



\Co\T 



> 4^\S 5 cos20|. 



(IV.46) 



200 MeV 

However this inequality requires a large value of |£g|, for example for T ~ lOMeV it requires that \Cg\ > 140. 



• case II: a; = fc + Mr- , ft = 1, positive energy, positive helicity neutrinos: 



10 -ie /TV 
<5 5 VMeV / 



M 
"eV 



^9 T 



2 MeV J 



(IV.47) 



Where we have neglected terms of higher order in Mr. Because M ~ leV and 100 MeV ^> T ^> 1 MeV a 
resonance would only be available for /c ~ 10~ 16 MeV which is not a relevant range of momenta for neutrinos in 
the plasma. Therefore, positive helicity neutrinos mix with the vacuum mixing angle. 

• case III: lu = —k — , h = — 1, positive energy, negative helicity anti-neutrinos: 



Ah 
SM 2 



10 -i6 /TV 
5 5 VMeV/ 



M 
"eV 



£9 



T 



VMeV/ 



(IV.48) 



Again in this expression we have neglected higher order terms in M^- . A conclusion similar to that of case II above 
holds in this case. No resonance is available for relevant values of neutrino momenta within the temperature 
range in which these results arc valid. For the cases II and III the ratio \Ah/5M 2 \ -C 1 for all relevant values 
of the neutrino momentum within the temperature range in which these results are valid. Therefore, negative 
helicity antineutrinos mix with the vacuum mixing angle, just as positive helicity neutrinos. 



• case IV: u> 



-k 



Ml 

2k 



1, positive energy, positive helicity anti-neutrinos: 



A h £ /0.1TV k_ 
SM 2 ~ 5 5 VMeV / T 



£9 + 



2T 
MeV 



k 
T 



(IV.49) 



Where we have neglected higher order terms in 4™- . The position of the resonance in this case is obtained from 
that in case I above by the replacement £9 — > —£9, namely for fixed temperature the resonance condition is 
fulfilled at the value of k given by 



(MeV) 2 
8T 



Co 



C% + 16 65 cos 29 



200 MeV V 
T ) 



(IV. 50) 
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Again in the temperature range 1 MeV <C T <C 100 MeV there is a resonance if 65 cos 28 > (assuming that 
\S 5 cos26»| - 1). 

Just as in case I, if \Cg\ < 1 the non-local term dominates over the asymmetry contribution for \/Tk > lMeV 
and the resonance occurs for k ~ 25^/5$ cos 28 (MeV) 3 /T 2 . 

Cases III and IV reveal an interesting feature: only the asymmetry contribution changes sign between neutrinos 
and antineutrinos whereas the non-local (in space-time) term remains the same. 

Together these expressions confirm that if the lepton and neutrino asymmetries are of the same order as the 
baryon asymmetry namely 0.2 < |£g| < 0.7, then the non-local terms from the exchange diagrams dominate the 
self-energy for T > 3 — 5 MeV unless the neutrino in the plasma has a momentum k such that VkT <C 0.5 MeV. 

In summary, for m e <C T <C resonances occur in cases I and IV when h A < (A is the sign of the energy 
eigenvalue). For h A > (cases II and III) no resonance is available for neutrino momenta k ~ T ~ few MeV and 
mixing angle in the medium coincides with the vacuum value. 



2. m e ,m M < T < M w 



We use m M » 106 MeV and find the following simple expressions 



• case I: w = k ■ 



M 
2k 



h = — 1, positive energy, negative helicity neutrinos: 

4 



A ft 
SM 2 



0.4 x 10 12 / T y k 
I 5 \GeVj T 



4.83 - + 10- 3 £ 9 



(IV.51) 



In this case no resonance is available but for neutrinos with extremely low energy and not relevant for the 
plasma. For example for T ~ GeV only neutrinos with energy of a few eV would be potentially resonant, but 
this momentum range is not a relevant one for neutrinos in the plasma. For neutrinos with energy larger than 
a few keV the mixing angle effectively vanishes. Therefore, we conclude that in this temperature regime the 
mixing angle in the medium for negative helicity neutrinos vanishes. 



• case II: u> = k + 4^- , h = 1, positive energy, positive helicity neutrinos: 



SM 2 



10- J ( T \ 2 T 



M 



4.83 



T 



10^£c 



(IV.52) 



It is clear that for the relevant regime of neutrino momenta in the plasma |A^/<5Af 2 | <C 1. Hence the mixing 
angle in the medium coincides with the vacuum mixing angle. Thus the conclusion in this case is similar to that 
in the case described by eq. (|IV.47fl . namely positive helicity neutrinos undergo oscillations in the medium with 
the vacuum mixing angle. 



• case III: lu = —k — , h — —1, positive energy, negative helicity anti-neutrinos: 

2 m ,— n 2 - 



Aft 
SM 2 



10- 



T \ T 
GeVJ ~k 



M 



4.83 - + 10- 3 £ 9 



(IV.53) 



The result in this case is similar to that of case II above, negative helicity antineutrinos oscillate in the medium 
with the vacuum mixing angle. 



• case IV: u> = — k — 4n- 

2k 



h = 1, positive energy, positive helicity anti-neutrinos: 

4 



SM 2 



0.4 x 10 12 /TV k_ 
I 5 VGeVy T 



4.83 



T 



10- 3 £ £ 



(IV.54) 



The conclusion in this case is similar to that of the case described by eq. ljiV.51(l above, the mixing angle 
effectively vanishes and oscillations of positive helicity antineutrinos are suppressed in the medium in this 
temperature range. 
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Taken together the above analysis reveals that there is a resonance in the oscillation of negative helicity neutrinos 
and positive helicity antincutrinos (that is ft. A < 0) in the temperature range m e «T« m„ with a typical neutrino 
momentum fc ~ T ~ fewMeV. For m e , m (1 «T« Myy the mixing angle for negative helicity neutrinos and positive 
helicity antincutrinos (that is h A < 0) effectively vanishes in the medium, and in both temperature ranges positive 
helicity neutrinos and negative helicity antineutrinos undergo oscillations in the medium with the vacuum mixing 
angle. We cannot yet conclude that positive helicity neutrinos and negative helicity antincutrinos are sterile, before 
studying the corrections to the oscillation frequencies. 



C. Oscillation frequencies and time scales 

The oscillation time scale in the medium is given by 



r# J (M) 



1 



1 



uj{ h) (fc, A) - J 2 h) (k, A) Ei (fc) - E 2 (fc) + Su[ h) (k, A) - 5J 2 h) (k, A) 



(IV.55) 



where Ei 2 (k) = yk 2 + M 2 2 and A = +1 and A = —1 correspond to neutrino and antineutrinos respectively. 
The vacuum oscillation time scale is 



T v (k) 



1 



(IV.56) 



\E 1 (k)~E 2 (k)\ ' 

therefore, in order to understand the loop corrections to the oscillations time scales, it is convenient to study the ratio 



T v {k) 



r£\k,X) 



Suj{ h) (k,X)-dLo { 2 n >(k 1 X) 



1 



E x {k)-E 2 (k) 



Typical neutrino momenta in the plasma are ultrarelativistic, hence we approximate 

SM 2 5 5 lQ- n eV 



E 1 {k)~E 2 (k) 



furthermore to leading order in Gf we replace ui^ '(fc, A) rs Afc in the arguments of A e>IJt {ui). 

The term Su^ (fc, A) — Slj^ (^j ^) represents the correction to the vacuum oscillation time scale due to the medium 
effect. While the general form of these corrections are cumbersome, we can extract simplified expressions in three 
relevant limits. 

• I Resonant case: A, '^^' 2 ' fc ' > cos 29: 



2k 



2 (fc/MeV) 



(IV. 57) 



(IV., 



In section TlV Bl above we found that resonant flavor oscillations can occur only for A = +l,h = —1 and 
A = — l,h = +1. In both these cases we obtain, 



Suj[ h) (k,X)-SLu y 2 l) (k,X) 



SM 2 
"8fc^ 

SM 2 



A e (k) + A^k) + B e (k) + B^k) 



M 



2fc 2 ,1 + ^IMl-sin^). 



(IV.59) 



• II vanishing mixing angle 



A h (A£i, 2 ,fc) 
SM 2 



> 1: 



In this limit, 9m ^ ~ and neutrino flavor mixing is suppressed. In section TlV Bl above, we found that this occurs 
only for A = +l,/i = —1 or A = —l,h = +1. Furthermore, Ah(XEi j2 , fc) is always positive definite in both 
temperature limits considered here m e «T« and m e , «T« Mw- 
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In this case we obtain 



5cj[ h \k,X)-8uj y 2 a> {k,X) 



SM 2 



A e (k) + A^k) + B e {k) + B^k) 



+ sign(5M 2 ) 



A e (k) - A^{k) - B e {k) + B^{k) 



• III vacuum mixing: 



A h (\E li2 ,k) 



8M* 



< 1: 



In this limit 6m ~ 0. In section llV Bl above we found that this case occurs for A = +1, h = +1 or A 
In both these cases we find 



54 h \k,X)-Sui h \k,X) = 

2 



M 
4k 2 



A e (k) + A^k) + B e (k) + B„(fc) 
• cos Jrt A e (k) - A^k) + B e (k) - B„(k) 



(IV.60) 



-l,h = -1. 



(IV.61) 



We now study these simplified expressions in the different regimes of temperature and for the different helicities 
components. The most relevant cosmological regime corresponds to momenta of the order of the temperature, hence 
we will focus on the regime in which the non-local (in space-time) contributions from the exchange diagrams dominate 
over the lepton- neutrino asymmetries. Taken together these simplifications allow us to study the relevant cosmological 
range of neutrino energies in a clear manner. 



m e «T« nip 



M 



case I: wi, 2 (fe, A) = A (k + -^-) ; A = +l,h = —1 and A — -1, h = +1: 



As observed in section (IIV Bjl the non-local terms become dominant for vTk > lMeV which is of course 

consistent with the ultrarelativistic limit M /2k 2 <C 1. Furthermore in the temperature range of interest in this 
study, the factors A e>M (fc) and B e ^(k) are of the order G F kT 4 /M^ ~ lCT 9 (T/GeV) 4 k < k. Therefore, near 
the resonance which occurs when T 2 k ~ 25 \J 5§ cos 26>MeV 3 , the expression l|IV.59f> simplifies to 



Su[ h \k,X)-SoJ 2 Jl) (k,X) 



The ratio of the oscillation time scales <jiV.57() becomes 

T v {k) ^ 
T^ ] (k,X) ~ 



f2 



SM 
'^2k 



(1 - sin 26») . 



sin2 01 < 1 



(IV.62) 



(IV.63) 



Therefore, for small vacuum mixing angle there is a considerable slow down of oscillations. Resonant flavor 
mixing in the medium occurs on longer time scales than in the vacuum. 

For large neutrino energy, well outside the resonance region for T 2 fc > 25 VI<55COs2(9|MeV 3 , eq. 11 V. 4411 
indicates that \Ah(XEi^, k)/5M 2 \ ^> 1. In this high energy regime we find that 



5w[ h) (k,X)-Suj^ ) (k,X) » sign(<5M 2 ) 

= sign(5M 2 ) 



A e (k) - A^k) - B e {k) + B^k) 



^ 28 7T 2 G F k T 4 
4571 M 2 , 



7.9xlO-eVsign(.M 2 )^ (-^J 



(IV.64) 
(IV.65) 
(IV.66) 
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therefore neutrino oscillations are suppressed by a vanishingly small mixing angle and the ratio of time scales 
(|1V.57|) becomes 



T v {k) 



HT 3 ( kT 2 



\5 5 \ VMeV-' 



(IV.67) 



A considerable speed-up of oscillations occurs for k T 2 > 100 MeV 3 since then Tm\k,X) <C r v (k). In this case, 
off-resonance flavor mixing is suppressed not only by a small mixing angle in the medium but also by a rapid 
decoherence and dephasing of the oscillations. 



• case II: idi^(k, X) = X (k 



2k 



■) ; A = +l,h= +1 and A = -l,h= -1: 



The results of section HVBI fsee eq. (|IV.47fl 1 indicate that in this case \Ah(XEi^, k)/SM 2 \ <C 1, corresponding 
to the mixing angle in the medium being the same as in the vacuum. As a result, 

-2 



M 
M 2 



cos 2 e 



cos 26* 



4 k 2 

-6.1 x 10~ 29 eV cos26> 



A e (k) - A^k) + B e (k) - B^(k) 

14 7T 2 Gf k T 4 
45 V2 M%r 



( k 



( T V 



VMeV/ VMeV/ 



and the ratio of time scales l|IV.57(l becomes 

T v (k) 



(fe,A) 



10- 17 cos 26» [TV 
6~ 5 VMeV ) 



(IV.68) 

(IV.69) 
(IV.70) 

(IV.71) 



Therefore, for positive helicity neutrinos and negative helicity antineutrinos medium oscillations are the same 
as vacuum oscillations both in the mixing angle as well as in the oscillation time scales. In this regime of 
temperature positive helicity neutrinos and antineutrinos are sterile in the sense that these do not interact with 
the medium. 



2. m e , m M < T <C M\ 



M 

case I: wi, 2 (fc,A) = A(fc+ -^f-) ; A = +1, h = -1 and A = —1, h = +1: 

This case describes negative helicity neutrinos and positive helicity antineutrinos. Eq. (|IV.51|) indicates that 
in this case, \Ah(XEi t 2 7 k)/SAl 2 \ 3> 1 corresponding to vanishing mixing angle in the medium. Therefore, we 
obtain 



SLo[ h) (k 1 X)-5J 2 h) {k,X) 



sign (6M 2 ) 
sign(5M 2 ) 



A e (k) - AJk) - B e {k) + BJk) 



2 Gf k T 2 



3 72 

9.6 x 10~ 6 eV sign(5M 2 



M w 
( k 



VMeV 



t y 

GeVJ 



and the ratio of time scales is given by 

T v (k) 



r£\k,X) 



10 12 / kT 



\S 5 \ VGeV z 



> 1 



(IV.72) 

(IV.73) 
(IV. 74) 

(IV.75) 



Hence there is a considerable speed-up in the oscillation time scale in the medium. Again, in this case oscillations 
are strongly suppressed not only by a vanishingly small mixing angle but also by the rapid dephasing in the 
medium. 
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case II: wi, 2 (fc,A) = A (k + -^f-) ; A = +1, h = +1 and A = -1, h = -1: 



This case describes positive helicity neutrinos negative helicity antineutrinos. Eq. I|IV.52(I shows that in this 
case | A/ l (A£^i ! 2, A;)/5M 2 | <C 1, the mixing angle in the medium is the same as in the vacuum. We find for this 
case 



Scj[ h) (k,X)-6u ( 2 ' l) (k,X) 



M 
Ik? 

r 2 



cos 2$ 



A e (k) - A^(k) + B e (k) - B^k) 



M~ 2 Gp k T 2 ( m M 
—— cos 26 I — ^~ 

4 k 2 

1.5 x l(T 19 eV cos26> 



3 y/2 \M W 
( k 



- 1 ( T \ 2 
VMeVy \GeVj 

where we have taken M ~ 1 eV. The ratio of time scales in this case is given by 



T v (k) 



Ah) 



(M) 



1 



3 x 1CT 8 cos20 / T 



GcV 



1 



(IV.76) 

(IV.77) 
(IV.78) 

(IV.79) 



Again in this temperature regime we find that positive helicity neutrinos and negative helicity antineutrinos 
are almost "sterile" in the sense that neither the mixing angle nor the oscillation time scales receive substantial 
loop corrections. Thus the combined analysis of mixing angle and propagation frequencies in the medium in 
the temperature regime under consideration indicates that in medium corrections for positive helicity neutrinos 
and negative helicity antineutrinos are very small. These degrees of freedom are effectively sterile in that their 
dynamics is (almost) the same as in the vacuum. 



LAPLACE TRANSFORM AND REAL-TIME EVOLUTION 



The main purpose to obtain the Dirac equation in real time is to study the oscillations of neutrinos in the medium 
as an initial value problem. As described in section[n]this is achieved by adiabatically switching the sources 77,77 from 
t = —00 and switching them off at t = 0. The adiabatic switching on of the sources induces an expectation value, 
which evolves in the absence of sources for t > 0, after the external source has been switched off. It is convenient to 
write the effective Dirac eq. (|II.17|) in terms of spatial Fourier transforms. Using the results of the appendix (|A 3|) we 
find 



{ n °§t ~ ^ ' y Sab " Mab + E ^ dL 



Mk, t) + / alt 1 E Qb (fc, t — t')L Mk, = ~Va(k, t), 



where the results of appendix IjA 3() yield 

dko 



Y.(k,t-t')=i I — ImE(fc,fc ) e- ifeo( *-*') , S(fe, fc ) = Z w (k, k ) + E z (k, k ) 

7T 



(V.l) 



(V.2) 



For an external Grassmann valued source adiabatically switched on at t 

Va(k,t) =Va(k,0) e e * 6{-t) , e- 



(V.3) 



— 00 and off at t = 
0+ . 

It is straightforward to confirm that the solution of the Dirac eq. (|V.1|I for t < is given by 

MKt<0) = Mk,0)e et . (V.4) 

Inserting this ansatze into the equation l|V.l|) it is straightforward to check that it is indeed a solution with a linear 
relation between ip a (k,0) and rj a (k,0). This relation can be used to obtain ip a (k,0) from rj a (k,0), or alternatively 
for a given initial value of the field at t — to find the source n a (k,0) that prepares this initial value. For t > the 
source term vanishes, the non-local integral in ea. (|V.l|l can be split into an integral from t = —00 to t = plus an 
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integral from t = to t. In the first integral corresponding to t < we insert the solution ea. l|V.4|) and obtain the 
following equation valid for t > 



d 



at 



H u ^7-7-M 6 ab -M ab + Z% d L 



ip b {k, t) + / dt' £ ah (£, t-t')L ij b (k, t') 



+ 30 



dk ImS ab (fc, fc ) 



7T fc 

This equation can be solved by Laplace transform. Introduce the Laplace transforms 

'"+ 00 dk ImS(fe,fe ) 



4>(k,s) 



dt e~ st tp(k,i) , S(fc,s) 



eft e" st £(fc,t) 



(V.5) 



(V.6) 



JQ J -oo n fco — is 

where we have used eq. ((V.2(l to obtain the Laplace transform of the self-energy, which leads to the analyticity relation 
(see Eq. (jnLU), 

S(fc,s) =£(fc,w = is-ie) . (V.7) 
In terms of Laplace transforms the equation of motion becomes the following algebraic equation 

(i 7 °s - i ■ fc) S ab - M ab + £*f L + X ab (k, s) i] Mk, s) = 1 1 7 ° S ab + 1 [£ o6 (fc, s) - S a6 (fc, 0)] L j ^ 6 (fc, 0) . 

(V.8) 

Consistently with the expansion of the self-energy in frequency and momentum up to order uj/M\y, and using eci. (|V.7fl . 
we replace the expression in the bracket in eq. I|V.8(I by 



1 r 

is 



dH{k,cj) 



s=0 



j=0 



E'(jfe,0) 



S b(fc, s) - S Qi) (fc,0) 

Using the representation QUI. 15(1 for the real part of the self-energy, the eq. I(V.8(I can be written as 
(7 i s - 7 • kj 1 - M + 7°A(s) L - 7 • kB(fe) L xf b (k, s) = i 7 [1 + A'(0) L] i> b (k, 0) 



where 

~ d A 

A(«)=A(w = ia) ; A'(0) = — 

The real time evolution is obtained by the inverse Laplace transform, 

r ds — . 
lKM)= / it-- i>{k,s)e si , 



j=0 



(V.9) 



(V.10) 



CV.11) 



(V.12) 



where T is the Bromwich contour in the complex s plane running parallel to the imaginary axis to the right of all 
the singularities of the function ip(k, s) and closing on a large semicircle to the left. We now follow the same steps as 
in section llVl namely projecting onto right and left components and onto helicity eigenstates. After straightforward 
manipulations we arrive at the following set of equations 

-(s 2 + k 2 ) 1 + (is - h(k)k)(A + h(k)B) - M 2 j xji L {k, s)=ij°M ip R (k, 0) + i(is - h(k)k) D ip L (k, 0) , (V.13) 

where D = 1 + A'(0), and 



ipR(k,s) = - 



is + h(k)k 



(V.14) 



We now follow the same steps as above to separate the Dirac and flavor structures by introducing the flavor doublets 
ip(k,s), £(fc,s) which are the Laplace transform of the flavor doublets cp(k,t), t) introduced in the expansion of 
the Dirac spinors in eas. ((rV.12() - ((iV.13(l . projecting onto the Weyl spinors eigenstates of helicity, the above equations 
become 

-(s 2 + k 2 )l + (is - hk)(A + KB) - M 2 1 (? (,i) (fc, s) = -i M^ h) (k, 0) + i (is - hk) D p {h) (k, 0) , (V.15) 
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is + hk ' 
' s 2 + k 2 



-M^ h) (k,s)+i^(k, 0) 



The solution to eq. (|V.15(I is obtained as 

^ h) (k, s)=^ h) {k,s) 



iM £ (h) (k,0) + i(is-hk)B (p {h) (k, 0) 



where the propagator is given by 



S h {k,s) 



1 



cth(k, s) + P(k, s) cos 29 



00 



-/3 h (k,s)sm26. 



(I') 



a 2 h (k,s)~ Pl(k,s)\ -/3 h (k,s)sm29 ( £ ) a h (k, s) - h {k, s) cos2^ ) 
in which it will prove convenient to introduce the following quantities 



cth(k, s) 



co 2 -k 2 + ^ S h (iu, k) - i (M 2 + M 2 ) 



(V.16) 



(V.17) 



(V.18) 



(V.19) 



LJ—is—te 



Ph(k,s) 



2 6M2 



cos 26 



A h (uj,k) 
SM 2 



+ sin 2 29 



LU—is — ie 



(V.20) 



The inverse Laplace transform ea. l|V.12() can be done straightforwardly, the singularities of <p^ (k, s) in the complex 
s plane are determined by the singularities of the propagator S^ h '(k, s). Up to the order in weak interactions considered 

here, these singularities are isolated poles along the imaginary axis at the positions s = — i u>£ given by eq. 

(|iV.36(l -l |IV.39(l . As a relevant example of the real time evolution, let us consider that the initial state corresponds 
to a wave-packet of left-handed electron neutrinos of arbitrary helicity h, with no muon neutrinos. This could, for 
example, be the case relevant for nucleosynthesis in which a neutron beta decays at the initial time. In this case 



<^ ( ' l) (fc,0) = u^{k) 



a [h) {k,o) 







(V.21) 



and wc hnd 



^ h \k,s) = vf\k) 



i{is-hk)[l + A'M\ ( a h (k,s)+f3 h (k,s) cos(20^ 



a 2 (k, S )-(3 2 (k,s) 



-0 h (k,s) sin(2£ i) ) 



(V.22) 



^(k,s) = -vW(k) 
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a 2 (k,s)- I3 2 (k,s) 
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a h {k,s)+l3 h (k,s) cos(20#>) 
-0 h (k,s) sin(2^ } ) 



(V.23) 



In order to avoid cluttering the notation, we have not included the frequency argument in the mixing angle in the 
medium 9$ but such dependence should be understood throughout. 

The term [cth{k, s) — (3h(k, s)] _1 features poles at s = —iu)i(k, ±) and the term [oth(k, s) + @h(k, features poles 
at s = —iu)2(k, ±). 

We neglect terms of order Gp T 4 /M^ ~ {T/M w ) 4 <C 1 since in the regime in which the approximations are valid 
T <C Mw- The residues of these poles are respectively 2wi,2(fe, ±) and the inverse Laplace transform yield within 
these approximations, 
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2. A 
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(V.25) 



For economy of notation, we introduce the following shorthand 



C 



(h) 

a, A 



s 



(h) 

a, A 



26^{^\k,X)) 



(V.26) 



for a = 1,2, helicity components h = ± and positive and energy components A = ±. 

The above expressions yield a direct comparison with the usual oscillation formulae in the literature. To leading 
order we set uji^ik, A) ps Afc (1 + M 2 2 /2k 2 ) in the prefactors in the expressions above thereby neglecting terms of 
order Gf- Since the dependence of the mixing angle on the frequency and momentum appears at order Gf, we 



can set u}i^{k, A) « A k in the argument of the mixing angles, therefore to leading order C\ 
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With these approximations, for an initial left handed electron state of helicity h = =F wc find 
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(V.28) 



The exponentials e^ 1 "^)* correspond to positive energy (— ) neutrino and positive energy (+) antineutrino compo- 
nents respectively. Therefore, the expressions above reveal that for negative helicity the relevant components in the 
relativistic limit correspond to positive energy neutrinos, while for positive helicity they correspond to positive energy 
antineutrinos. 

The upper component of the expressions above correspond to wavepackets of negative and positive helicity respec- 
tively for a left handed electron neutrino, namely i/^(k,t) while the lower components correspond to a left handed 
muon neutrino, namely z/J(fc,t). 

For the right handed components £ the leading order can be obtained by setting M«M1 thereby neglecting terms 
of order 6M 2 /M . This approximation yields, 
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(V.29) 



The upper and lower components correspond to wavepackets for right handed negative and positive helicity electron 
and muon neutrinos respectively. 

From the expression (|V.27(I we can read off the probability for relativistic left handed, negative helicity electron 
and muon neutrinos as a function of time, 



\v- L (k\t)\ 2 = K(fc)ni-sin 2 [2 m (fc)] sin 2 ( 1 Aoj~(k,+)t) +0 



Mj 2 
k 2 



K L (k,t)\ 2 = \v-(k)\ 2 {sm 2 [2 6-(k)] sin 2 (- Ao;-(k,+)t) +0 



'Mi* 
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(V.30) 
(V.31) 
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The probability for left handed relativistic positive helicity electron and muon antineutrinos as a function of time are 
read off from eq. l|V.28fl 

I<l(M)| 2 = kt(fc)| 2 jl-sin 2 [2#+(--fc)] sin 2 QA^+(fc,-)^+0^^) j , (V.32) 
I<l(M)| 2 = k + (£)| 2 {sin 2 [2 0+(-fc)] sin 2 Q Ac+(fc,-)^+0^^| , (V.33) 

where 

Aw ± (fc,±) = ^- + 5ut(k,±)-8ut(k,±) . (V.34) 
2k 

The corrections S(j^(k, ±) had been studied in detail in section II V CI above. Finally, eq. I|V.29|) determines the 
probabilities of finding right handed positive and negative helicities neutrinos as a function of time. This equation 
makes manifest that this probability is suppressed by a factor M/k with respect to the left handed component, indeed 
it is the mass term that is responsible for generating a right handed component from a left handed one and must 
therefore be suppressed by one power of the ratio M/k. For a typical neutrino momentum k > 1 MeV this suppression 
factor is of order 10~ 6 . Eas. (|V.24(l - irV25l) provide a complete field theoretical description of oscillations in real time. 
E as. IjV. 30(1 - IjV. 31(1 are obviously reminiscent of the familiar oscillation equations obtained in the simplified quantum 
mechanical two level system, however there are some important aspects that must be highlighted, namely, the field 
theoretical formulation introduced here led directly to these oscillation formulae in terms of the mixing angles in 
the medium and the correct oscillation frequencies that include the quantum loop corrections. Furthermore, the 
oscillation formulae obtained above reveal the nature of the approximations and allow a consistent inclusion of higher 
order effects as well as describe the oscillation of all helicity components as well as the dynamics of the right handed 
component. The usual oscillation formula obtained within the single particle quantum mechanical description emerge 
cleanly in suitable limits and the nature of the corrections is readily manifest. 



VI. CONCLUSIONS, AND FURTHER QUESTIONS: 



We have provided a systematic treatment of neutrino oscillations and mixing directly from quantum field theory in 
real time in a regime of temperature and density relevant for early Universe cosmology prior to big bang nucleosynthe- 
sis. While we have focused on two flavors (electron and muon) of Dirac neutrinos the formulation can be generalized 
straightforwardly to more flavors and to Majorana-Dirac mass matrices. 

We have obtained the medium corrections to the dispersion relations and mixing angles of propagating neutrinos. 
Implementing methods from real time non-equilibrium quantum field theory at finite temperature and density we have 
systematically obtained the equations of motion for the neutrino field and studied the real time evolution as an initial 
value problem. The major advantage of this approach, as compared to the usual approach based on single particle 
quantum mechanics is that it consistently and systematically includes the medium corrections to the dispersion 
relations and mixing angles directly into the real time evolution and treats left and right handed fields and both 
helicity components on equal footing. We have argued that collisional relaxation yields thermalization of neutrinos in 
flavor eigenstates for temperature T > 5 — 10 Mev for which the relaxation time scale via weak interactions is shorter 
than the oscillation time scale. Assuming the validity of this argument, we obtained the neutrino self-energies up to 
one-loop including the asymmetries from leptons, neutrinos, hadrons and quarks, as well as non-local (in space-time) 
terms arising from the expansion of the self-energy loop in the external frequency and momentum. We have consider 
these non-local terms up to leading order in co/Myy ■ k/Mw since these terms are of the same order of or larger 
than the contribution from the asymmetries for T > 5 MeV if all asymmetries are of the same order as the baryon 
asymmetry. This is yet another indication that this is an important temperature regime in the early Universe. 

Our main results are summarized as follows: 

• Implementing the methods from non-equilibrium real-time field theory at finite temperature and density we 
obtained the equations of motion for the neutrino fields in linear response. This formulation includes consistently 
the self-energy loop corrections to the dispersion relations and mixing angles in the medium and treat left and 
right handed fields with both helicity components on equal footing. 

• We have focused on a temperature regime prior to nucleosynthesis T > 5 — 10 MeV in which we argued that 
neutrinos are thermalized as flavor eigenstates. We studied two different temperature regimes: m e <C T <C 
within which we have shown that there is the possibility of resonant oscillations of test neutrinos, and m e , m M <C 
T <C M\\r within which the mixing angle for active neutrinos effectively vanishes. 
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• An expansion of the self-energy in terms of the neutrino frequency and momentum is carried out to lowest order in 
ui/M\y ; k/My/ thus extracting the leading non-local (in space-time) contributions. We find a new contribution 
which cannot be identified with an effective potential. The mixing angles and propagation frequencies in the 
medium are found to be helicity dependent. It has been recently pointed out[46j that a neutrino helicity 
asymmetry could be a very important ingredient in successful see-saw models of thermal leptogenesis. 

• If the lepton and quark asymmetries are of the same order as the baryon asymmetry in the early Universe, 
we have shown that the non-local (in space-time) terms in the self-energies dominate over the asymmetry for 
typical energies of neutrinos in the plasma for T > 3 — 5 MeV. 

• The oscillation time scale in the medium is slowed-down near the resonance, becoming substantially longer than 
in the vacuum for small vacuum mixing angle. For high energy neutrinos off-resonance the mixing angle becomes 
vanishingly small and the oscillation time scale speeds-up as compared to the vacuum. At high temperature, 
in the region T 3> m ei m^ the mixing angle for active neutrinos effectively vanishes and there is a considerable 
speed-up of oscillations, which are then suppressed by a vanishingly small mixing angle and a rapid dephasing. 

• We have obtained the general equations of motion for initially prepared wave packets of neutrinos of arbitrary 
chirality and helicity. These equations reduce to the familiar oscillation formulae for ultrarelativistic negative 
helicity neutrinos, but with the bonus that they consistently include the mixing angles and the oscillation 
frequencies in the medium. These equations not only yield the familiar ones but also quantify the magnitude 
of the corrections. Furthermore these equations also describe the evolution of right handed neutrinos (of either 
helicity) which is a consequence of a non-trivial mass matrix and usually ignored in the usual formulation. 

Further questions: The next step in the program is to provide a solid assessment of the reliability of the 
argument that suggests that for T > 10 MeV neutrinos thermalize as flavor eigenstates. The phenomenon of slow- 
down of oscillations near a resonance found in this work, suggests that this argument is consistent at least near 
resonances. We are currently studying the relaxational dynamics by including two loop corrections to the self-energy 
which account for collisional processes in the medium. We are also studying how to obtain the kinetic equations that 
describe simultaneously oscillations and relaxation systematically within a field theory approach that accounts for the 
subtle aspects of the Fock representation necessary to understand the concept of the flavor distribution functions [37I. 
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APPENDIX A: REAL-TIME PROPAGATORS AND SELF-ENERGIES. 

1. Fermions 

Consider a generic fermion field f(x,t) of mass m/. The Wightmann and Green's functions at finite temperature 
are given as 

i S> p (2-5>, t-t') = (f a (x, t)-f p {x", t')) = I £ e^-*') iS>f,& t - t'), (A.l) 

p 

iS<^(x-x',t-t') = -(7 / ,(/ 1 f)/ a (? 1 t)) = ^e^i^(f,t-t') ) (A.2) 

p 

where a, [3 are Dirac indices and V is the quantization volume. 

The real-time Green's functions along the forward (+) and backward (— ) time branches are given in terms of these 
Wightmann functions as 

(/£ +) 0?.*)7j, + V ,0) = i S ++ (x - x 1 ,t - t') = i S > (x - X 1 ,t - t') Q(t-t') + i S<(x- x',t-t') 9(t'-t)(A.3) 
(fi +) (x,t)J ( p ~\x',t')) = iS + -{x-x\t-t')=iS < {x-x\t-t') . (A.4) 
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At finite temperature T, it is straightforward to obtain these correlation functions by expanding the free fermion fields 
in terms of Fock creation and annihilation operators and massive spinors. In a CP asymmetric medium, the chemical 
potential pj for the fermion / is non-zero. Particles and anti-particles obey the following Fermi-Dirac distribution 
functions respectively 

N fW= efro-M/U + i ' ^ (po)= e (po+J)/T + 1 - ( A - 5 ) 
The fermionic propagators are conveniently written in a dispersive form 

/oo 

dpo (Z,f,&Po) e~ ipo (*"*') , (A.6) 

/OO 

dpoP< l3 (p,Po)e- i ^ t - t '\ (A.7) 
-oo 

where we have 



7°Po - 7 -p + rnf . 7° p - 7 • p '+ m/ 



P>>>o) = ' 2pQ " J [1 - N f {p )] S( P0 Wp ) + ! ^ ^ J A^-Po) <Kpo + w p ) , (A.8) 

P<>>o) = 7 ° P ° "J ' P " + ^ ^/(Po) *(Po - <* P ) + 7 " P ° ~J ' ^ + m/ [1 - fo(-po)] gfa, + , (A.9) 
zpo 2p 

with oj p = \J\p\ 2 + rrij. Using the relation Nf(—po) = 1 — Nf(po), we can write 

PZjftPo) = [1 - N f (p )} pl !0 (p,Po) , (A.10) 

pZ,/3&Po) = Nf( Po ) pIjMpo) , (A.11) 
where the free fermionic spectral density p^(p,Po) is given by 

P J (p,Po) = t; — ^(Po-^ P ) + T — (5(po + w p ), (A.12) 

P± = 7° =F 7 ' P i m / • (A. 13) 

2. Vector Bosons 

Consider a generic real vector boson field A^ix, t) of mass M. In unitary gauge, it can be expanded in terms of 
Fock creation and annihilation operators of physical states with three polarizations as 

A^x,t) = ^^^^ = Ja^ x e- iW ^e i ^ + al x e iW ^e- i ^ ; fc^, A (£) = 0, (A.14) 



fe 



where Wk — \J\k\ 2 + M 2 and fc^ is on-shell = (Wk , fc) . The three polarization vectors are such that 

Y j e»{k)el{k) = P^{k) = -U»-— r ) . (A.15) 



It is now straightforward to compute the Wightmann functions of the vector bosons in a state in which the physical 
degrees of freedom are in thermal equilibrium at temperature T. These are given by 

(A^{x,t)A v (x',t')) = iG>„(x-x , ,t-t'), (A.16) 

(A,(x',t')A^x,t)) = iG< v {x-x',t-t'), (A.17) 

where G K> can be conveniently written as spectral integrals in the form 

1 - f°° 

iGlJx-x',t-t') = YVM*-') / 4 e -Mt-t') [1 + N b (k )]p^(k ,k), (A.18) 

* - J — oo 

k 

iG< I/ (x-x',t-t') = VJ e *(---') / dioe-^C-''' N b (k ) p^(k ,k) , (A.19) 
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where 



N b (k ) 



3 fco/T _ I 



and the spectral density is given by 



2W k 



P»u(k) S(k - W k ) - P^i-k) 5(k + W k ) 



(A.20) 



(A.21) 



In terms of these Wightmann functions, the real-time correlation functions along the forward and backward time 
branches are given by 

= iGl v (x-x',t-t')Q(t-t')+iG< u (x-x',t-t')Q(t'-t), (A.22) 
(4+>(f,<)4-)(f',0) = iG<„(x-x',t-t'). (A.23) 

For the charged vector bosons, the correlation functions can be found simply from those of the real vector boson fields 
described above by writing the charged fields as linear combinations of two real fields A 1,2 , namely 



W±{2,t) = ^=[Al(x,t)±iAl(x,t)] . 

It is straightforward to find the correlation function 

(W+(x,t)W-(x',t')) = G> v (x-x' ) t-tf) , 
and similarly for the other necessary Wightmann and Green's functions. 



(A.24) 



(A.25) 



3. Retarded Self-Energies for Charged and Neutral Current Interactions 

The diagrams for the one-loop retarded self-energy from charged current interactions are displayed in Fig. A 
straightforward calculation yields for the charged current contribution the following result 

t-t') = l JLRY [i S ++ (x -x',t- t') i G++(x -x',t- t') -iS < {x-x',t- t') i G< v {x ~x',t~ t')] Y L 

(A.26) 



with 



R = 



1+7 5 



L = 



1 -7 5 



A similar result is obtained for the neutral current contribution to the self-energy by simply replacing g/V% — > 
g/2cos9 w and M\y — > Mz = M\y/ cos9 w . 

Using the representation of the fermion and vector boson propagators given above the retarded self-energy l|A.26|) 
can be written as 



5jret (x X , t t ) 



I 

V 



E 



dkn R 



Z w (k, k ) + Zz(k, ho)] L jHs-*') e -«o(t-f ) ( t _ t ' )j ( A _27) 



The contributions from charged and neutral vector bosons are given by 



^ 2 ' d 3 q 



~*h (/•-./.•.,) = / / dp I dq 6(po + qo-ko) J» ?,Po) P™{k)i>l </" > ~" ') --V.,-(/>n'l + X-,(«/u)] . (A.2N) 



Sz(fc,fc ) 



d 3 q 

4 cos 2 9 W J (2tt) 3 
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d Po / dq Q 5{ Po +q -k Q ) ^ pf{k-q,p Q ) p*(q,q ) Y [1- N f (p )+N b {q )] , (A.29) 



where pw,z(q,qo) are the vector boson spectral densities given by ea. l)A.21(l with M = Mw,z respectively. It is clear 
that Y*w,z(k, ko) corresponds to a vector-like theory. 
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x't') 




FIG. 2: Retarded self-energy for charged current interactions. The wiggly line is a charged vector boson and the dashed line 
a lepton. The labels (±) correspond to the forward (+) and backward ( — ) time branches. The corresponding propagators are 
i S ' (x — x , t ~ t') and i (x — x', t — t') for leptons and charged bosons respectively. 



Using the integral representation of the function @{t — t'), the retarded self-energy can be written in the following 
simple dispersive form 



E ret (x-x , ,t-t') = ^-J2 I T- e iHs ~ £,) e-W-V R Z w (k,u) + X z (k,Lu) 



V 



du> 



2tt 



(A.30) 



^w,z 



(k,k ) 



fco — lu — ie 



(AM) 



Sw,z(fe)W) = / dko 
where e — > + . Hence, from the above expression, we identify 

^wz{k,Lu) — — Im'Ew z{k,uj) ■ (A. 32) 

TT 

Furthermore, since R(±rrif)L — 0, the factor ±mj in the free fermionic spectral density defined in Eqns. ljA.12f) 
and ljA~13jl can be ignored when we compute R Y,w{k, w) + lu) L. The signature of the the fermion mass m/ 

is only reflected in the factors of w p in the spectral density. Ignoring the factor rrif from now on, the fermionic spectral 
density is proportional to the 7 matrices only and does not feature the identity matrix or 7 . Therefore, there is the 
following simplification 



R Y,w{k,oj) + Y,z{k,uj) L= u>) + ^z{k, lu) L 



(A.33) 
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